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Abstract: We calculate the relic abundance of thermally produced neutralino cold dark 

matter in the general 19 parameter supergravity (SUGRA-19) model. A scan over GUT 

scale parameters reveals that models with a bino-like neutralino typically give rise to a dark 

matter density Q^oh 2 ~ 1 — 1000, i.e. between 1 and 4 orders of magnitude higher than 

Al 

the measured value. Models with higgsino or wino cold dark matter can yield the correct 
relic density, but mainly for neutralino masses around 700-1300 GeV. Models with mixed 



bino-wino or bino-higgsino CDM, or models with dominant co- annihilation or ^4-resonance 
annihilation can yield the correct abundance, but such cases are extremely hard to generate 
using a general scan over GUT scale parameters; this is indicative of high fine-tuning of 
the relic abundance in these cases. Requiring that m^o ~ 500 GeV (as a rough naturalness 
requirement) gives rise to a minimal probably dip in parameter space at the measured CDM 
abundance. For comparison, we also scan over mSUGRA space with four free parameters. 
Finally, we investigate the Peccei-Quinn augmented MSSM with mixed axion/axino cold 
dark matter. In this case, the relic abundance agrees more naturally with the measured 
value. In light of our cumulative results, we conclude that future axion searches should 
probe much more broadly in axion mass, and deeper into the axion coupling. 
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1. Introduction 

The cosmic abundance of cold dark matter (CDM) has been recently extracted at high 
precision by the WMAP collaboration Q. An analysis of WMAP seven year data implies 

Q CDM h 2 = 0.1123 ± 0.0035 (68% CL), (1.1) 

where Q = p/p c is the dark matter density relative to the closure density, and h is the 
scaled Hubble constant. Since no particle contained within the Standard model of particle 
physics has exactly the right properties to constitute CDM, new matter states from physics 
beyond the SM are required for explanation ||. 

The most popular explanation of CDM comes from models including weak scale super- 
symmetry (SUSY)||. SUSY models are invoked to cure the hierarchy problem contained 
within the SM. SUSY models also receive some indirect support from experiment H in that 
the gauge couplings, measured at the scale Q = Mz, meet nearly at a point when run to 
Mgut via the renormalization group equations (RGEs) of the Minimal Supersymmetric 
Standard Model, or MSSMQ. Within i?-parity conserving SUSY theories, there are sev- 
eral dark matter candidates, including the lightest neutralino Xv ^ ne gravitino G and the 
axion/axino admixture aa from the Peccei-Quinn||, |7| augmented MSSM (or PQMSSM).|8| 

Most analyses focus upon the lightest neutralino Xi as a CDM candidate, since it can 
be classified as a weakly interacting massive particle, or WIMP. WIMP particles such as 
X? are especially compelling due to the so-called "WIMP miracle" . It is assumed the DM 
particle is present in thermal equilibrium at some early time in the universe's history, and 
then the relic abundance can be found by solving the Boltzmann equation for its number 
density. Particles with weak scale mass, which interact with weak scale strength, should 
give a relic abundance roughly in the vicinity of what is measured. 

However, detailed calculations of the neutralino relic abundance in models such as 
mSUGRA || (or CMSSM) find that throughout parameter space, the neutralino abun- 
dance is usually much higher than the measured value, and typically by 1-2 orders of 
magnitude [[1C|]. Only in very narrow regions of parameter space is the measured CDM 
abundance found: the stau co-annihilation region[12], the hyperbolic branch/ focus point 



region with mixed higgsino DM [13], the A-resonance region at large tan/3[|14|], and perhaps 
at the light Higgs resonance[15|. 

That the mSUGRA parameter regions where the predicted abundance agrees with 
measurement are quite narrow is indicative of a needed fine-tuning of fundamental model 
parameters to achieve the required value. The fine-tuning is a measure of the slope of the 
surface of Q^oh 2 as a function of model parameters aj. 

Quantitatively, the fine-tuning with respect to variation in the parameter aj is defined 
as 

s abgnw^ (L2) 

dlogat 
The overall fine-tuning can be obtained by combining the A ai in quadrature. The degree 



of fine-tuning of the relic abundance in mSUGRA has been calculated in Ref's [16] an d 



1 17]. Both groups find that the regions of agreement between theory and experiment are 



highly fine-tuned. Ref. |17J claims also that the very large tan/3 (~ 50) region is not 
very fine-tuned, since the pseudoscalar Higgs width Ta becomes very large (owing to the 
large 6-quark Yukawa coupling) , and efficient neutralino annihilation can occur all over the 
itiq vs. m 1 / 2 plane. 

In Ref. |18| , we proposed mixed axion/axino cold dark matter in the Peccei-Quinn 
augmented mSUGRA model. By requiring a re-heat temperature in the 10 6 — 10 8 GeV 
range, a MeV-scale axino mass and a large value of f a ~ 10 11 GeV is favored. As opposed 
to the high fine-tuning which is present in mSUGRA for neutralino CDM, if instead mixed 



axion/axino CDM is invoked, with a ~ MeV scale axino IS], then the fine-tuning needed to 
achieve a theory-experiment match is relatively low. The lowest fine-tuning is observed for 
the case where there is a nearly equal mix of axion and thermally produced axino DM0. 

A major criticism of the result found in Ref. [|l7]] is that the fine-tuning analysis of 
neutralino CDM is restricted to the mSUGRA model. It might be possible that if one 
opens up the parameter space to include models with non-universality, then many more 
possibilities may arise to gain a theory-experiment match in the dark matter relic density. 

In this paper, we calculate the relic abundance of neutralino CDM in a 19 parameter 
version of the MSSM: the SUGRA-19 model. Another recent calculation of the neutralino 
CDM abundance in a 19 parameter MSSM has been performed by Berger et a/jpof]. 1 In 
that analysis, the so-called pMSSM model is examined, which is a version of the MSSM 
with minimal CP and flavor violating terms. The 19 soft SUSY breaking and other pa- 
rameters are input at the weak scale, and it is required that Q^oh 2 < 0.11. In the analysis 
presented here, we will also adopt a 19 parameter MSSM, but will insist on adopting soft 
SUSY breaking parameters at the GUT scale. The reason is that the unification of gauge 
couplings provides circumstantial evidence that the MSSM, or MSSM plus gauge singlets 
(or additional SU(5) multiplets) is the correct effective field theory between M wea k and 
Mqut- The pattern of soft term running as provided by the MSSM RGEs gives us some 
guidance as to how likely are various sets of weak scale parameters. For instance, by sam- 
pling over weak scale parameters including the superpotential fi term, one easily generates 
low \x solutions which lead to higgsino-like neutralinos. However, sampling over GUT scale 
parameters, it is not easy to generate low /x solutions with higgsino-like dark matter. The 
reason is that (at tree level) ^ 2 {m wea k) — —mjj (M wea k), and it is rare that m 2 H runs 
to just barely negative values, which would lead to a small [i parameter. In addition, in 
distinction with Ref. p0(| , we W1U n °t restrict ourselves to requiring Q^oh 2 < 0.11, but will 
instead allow all values of thermal neutralino abundance. 

In Sec. ^, we provide an overview of our scanning methodology and parameter space 
values, and our calculation of the neutralino relic abundance. In Sec. |3|, we present 
results of the neutralino relic density from a scan over the SUGRA-19 model. We find 
that solutions with bino-like neutralinos give a large overabundance of CDM. Models with 
higgsino-like or wino-like neutralinos tend to give too little CDM, unless m^o ~ 1 TeV, 
with observable sparticles likely beyond LHC reach. The value O^-o/i 2 ~ 0.11 lies at the 
minimum probability between these two extremes. Since our scan depends on the bias as 

1 Similar analyses using weak scale SUSY parameters and requiring a thermal abundance of neutralino 
CDM in accord with WMAP are given in Ref's |I], ||, ||. 



to how we sample the GUT scale parameters, we present results for both a linear as well as 
a log scan; the latter scan accentuates small soft SUSY breaking parameter values. In Sec. 



3.3, we present for comparison a similar analysis done for the mSUGRA model. In Sec. |4|, 
we work instead within the Peccei-Quinn augmented MSSM, or PQMSSM, which leads to 
mixed axion/axino CDM if the axino is assumed to be the lightest SUSY particle (LSP). 
The PQMSSM requires four additional parameters: the PQ breaking scale f a and initial 
mis-alignment angle 9i, the axino mass ma, and the re-heat temperature of the universe 
Tr. In this case, the PQMSSM with 23 free parameters can more easily yield the measured 
relic abundance. We evaluate the favored range of the additional PQMSSM parameters. 
Our conclusions are presented in Sec. pi 

2. Relic density in the SUGRA-19 model 

For our calculations, we adopt the Isajet 7.80[24, 25 1 SUSY spectrum generator Isasugra. 



Isasugra begins the calculation of the sparticle mass spectrum with input DR gauge cou- 
plings and fb, f T Yukawa couplings at the scale Q = Mz (ft running begins at Q = m t ) 
and evolves the 6 couplings up in energy to scale Q = Mqut (defined as the value Q where 
Q\ = 92) using two-loop RGEs. At Q = Mqut, the SSB boundary conditions are input, 



and the set of 26 coupled two- loop MSSM RGEs [26] are evolved back down in scale to 
Q = Mz- Full two- loop MSSM RGEs are used for soft term evolution, while the gauge 
and Yukawa coupling evolution includes threshold effects in the one-loop beta-functions, 
so the gauge and Yukawa couplings transition smooothly from the MSSM to SM effective 
theories as different mass thresholds are passed. In Isajet 7.80, the values of SSB terms 
which mix are frozen out at the scale Q = Msusy = J m i L m t R 1 while non- mixing SSB 
terms are frozen out at their own mass scale f25|]. The scalar potential is minimized us- 
ing the RG-improved one-loop MSSM effective potential evaluated at an optimized scale 
Q = Msusy which accounts for leading two-loop effects [27]. Once the tree-level sparticle 



mass spectrum is computed, full one-loop radiative corrections are caculated for all sparti- 
cle and higgs boson masses, including complete one- loop weak scale threshold corrections 
for the top, bottom and tau masses at scale Q = Msusy- These fermion self-energy terms 
are critical to evaluating whether or not Yukawa couplings do indeed unify. Since the GUT 
scale Yukawa couplings are modified by the threshold corrections, the Isajet RGE solution 
must be imposed iteratively with successive up-down running until a convergent sparticle 
mass solution is found. For most of parameter space, there is excellent agreement between 
Isajet and the SoftSUSY, SuSpect and Spheno codes, although at the edges of parameter 



space agreement between the four codes typically diminishes [28]. We adopt the Isasugra 
non-universal SUGRA parameter space. 

We will implement at first a linear scan over the following parameters. 

• Gaugino masses: M 1 , M 2 , M 3 : - 3.5 TeV 

• First/second generation scalar masses: mQ l , mu X i rn D 1 - I tul^ tuex'- — 3.5 TeV, 

• Third generation scalar masses: rriQ 3 , mu 3 , mD 3 , thl 3 , rriE 3 '- — 3.5 TeV, 



• Higgs soft masses: m# u , mjf d : — 3.5 TeV, 

• trilinear soft terms: At, Aj,, A T \— 3.5 TeV — > 3.5 TeV, 

• ratio of weak scale Higgs vevs tan (3 : 2 — 60. 

We adopt a common mass for first and second generation scalars so as to avoid SUSY 
FCNC processes. 

To gain an acceptable sparticle mass solution, we will require: 

1. the lightest SUSY particle (LSP) is the neutralino x?> 

2. the lightest chargino, if non-wino-like, obeys the LEP2 limit vfiy x > 103.5 GeV, 

3. the lightest chargino, if wino-like, obeys the LEP2 limit m-^ > 91.9 GeV, 

4. the light Higgs mass obeys the LEP2 limit rrih > 111 GeV (where we allow a roughly 
3 GeV uncertainty in the theory calculation as applied to the actual limit where 
m h > 114.4 GeV. 

For each acceptable solution, we calculate the neutralino relic density Q^oh 2 using 
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the IsaReD|29|] program. IsaReD calculates all relevant neutralino annihilation and co- 
annihilation reactions, as obtained using CalcHEP, and then calculates the relativistic 
thermally-averaged (co)-annihilation cross sections times relative velocity. Once the freeze- 
out temperature is determined, then the relic density at the present time is found by 
integrating the Boltzmann equation as formulated for a FRW universe. 

Since we assume the neutralino to be in thermal equilibrium, our relic density results 
do not explicitly depend on the value of the re-heat temperature of the Universe Tr after 
inflation. However, we must assume Tr > Tj ~ m^-o/20 so that Tr is above the neutralino 

freeze-out temperature. Further, if Tr ~ 10 10 GeV, then thermal production of gravitinos 
in the early Universe, followed by decays to the LSP, will overproduce neutralino dark 
matter. Hence, for our neutralino CDM relic density calculations, we must assume here 
that 

• m^o/20 ~ Tr ~ 10 10 GeV. 

3. Results for neutralino cold dark matter 

3.1 Linear scan over SUGRA-19 parameters 

Our first results from a linear scan over the above SUGRA-19 parameter space is shown 
in Fig. [y, in the Q^oh 2 vs. m^o plane. 2 The various solutions are color coded according 
to the gaugino/higgsino content of the neutralino. In the notation of Ref. ||, if the 
bino-component \v\ | > 0.9, then the neutralino is labeled as bino-like (blue diamonds); 
if the wino-component \v% \ > 0.9, then it is labeled wino-like (purple x); if the higgsino 



2 A qualitatively similar plot was generated by Gelmini et aLfBOl (Fig. 7) using an MSSM model with 9 
input weak scale parameters. 



components yv\ + v 2 > 0.9, then it is labeled as higgsino-like (red squares). If the 
neutralino falls into none of these categories, then it is labeled as "mixed" DM: (orange 
circles). 

We see from Fig. |l| that the bino-like neutralinos tend to populate the region with 
Q^oh 2 S> 0.1, i.e. usually about 2-3 orders of magnitude too high. For low values of m^o, 

A 1 A 1 

the abundance tends to be more like 3-5 orders of magnitude too high. A few bino-like 
points do tend to make it into the Q^oh 2 ~ 0.1 region; these solutions tend to come form 
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various co- annihilation or resonance annihilation processes. To obtain the required relic 
abundance via co-annihilation, the LSP-NLSP mass gap must be tuned to just the right 
value. To obtain the required relic abundance via resonance annihilation, the LSP mass 
must be adjusted to be close to half the mass of the resonance. These co-annihilation and 
resonance annihilation points are quite hard, but not impossible, to generate using our 
random scan over GUT scale parameters. 

The higgsino-like and wino-like CDM bands also show up as distinct lines, typically 
with ST-o h? too low by 1-2 orders of magnitude unless m^o ~ 800 — 1200 GeV. The wino-like 
band is relatively well-populated, as this just requires Mi to be the lightest of the gaugino 
masses at the weak scale. The higgsino-like band is relatively less populated, showing 
that higgsino-like CDM is rather fine-tuned if one starts with GUT scale parameters, as 
mentioned in Sec. [j]. The points with the lowest population are those with mixed bino- 
higgsino-wino CDM. These "well-tempered neutralino |3l]]" points most naturally tend to 
populate the f^o/1 2 ~ 0.1 line, but they do require a fine-tuning to avoid a bino, wino 
or higgsino dominance. Especially at low m^o, relatively few solutions are found with 
O^o/i 2 ~0.1. 

To apprehend more clearly the dark matter probability distribution after our linear 
scan of SUGRA-19 parameter space, we project the model points listed in Fig. Ill as a 
histogram onto the O^-o/i 2 axis in Fig. |2|a). Here we see the most probable value of Qj,i no h 2 
is ~ 10 — 100 for bino-like dark matter (blue histogram), while the most probable value 
for wino-like dark matter is Q W i no h 2 ~ 0.005 — 0.05. The dip between these two cases 
is partially filled in by cases of bino, higgsino or wino, or a mixture, with the minimum 
probability lying around fi^o/1 2 ~ 0.2 — 0.4, i.e. just above the measured value. 

A large number of the wino and higgsino dark matter solutions with O^o/i 2 ~ 0.1 
come from models with very heavy neutralinos: m^a ~ 800 GeV. If the x? is the LSP, as 
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assumed here, then all other sparticles are heavier- and usually much heavier- than this 
value, and will likely lead to solutions with high electroweak fine-tuning [32]. In Fig. |2|&)., 
we plot the number of model solutions from the SUGRA-19 scan versus m^o, where in 
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addition we require (somewhat arbitrarily) that m^o < 500 GeV, so the solutions are not 
too fine-tuned with regard to electroweak symmetry breaking. In this case, the higgsino 
and wino LSP models which naturally give Jl^o/1 2 ~ 0.1 are all excluded. The maxima of 

Al 

bino-like solutions has moved up slightly to fibino^ 2 ~ 50 — 100, while wino-like solutions 
peak at Q w i no h 2 ~ 0.01. The minimum of the probability distribution lies very close to 
the measured value Ocdm^ 2 ~ 0.1. With sparticle masses generally at the TeV or below 
scale, the measured relic density lies at the least likely value as predicted by the SUGRA-19 
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Figure 1: Thermal abundance of neutralino cold dark matter from a linear scan over the SUGRA- 
19 parameter space. We plot versus the neutralino mass. Models with mainly bino, wino, higgsino 
or a mixture are indicated by the various color and symbol choices. There are 5252 points in the 
figure. 

model. In this case, it would be extremely fortuitous if the lightest neutralino of SUGRA 
theories was in fact the dark matter particle. 



3.2 Log scan over SUGRA-19 parameters 

We have interpreted our linear scan over SUGRA-19 parameter space in terms of a proba- 
bility distribution as to likely values of Q^oh? which would be obtained in super symmetric 
models. These results do depend on how we sample our GUT-scale parameter space. While 
it is impossible to know what the correct measure is for GUT-scale SUSY parameters, at 
least we can compare our linear scan against results using a different sampling measure for 
our parameter space. 

In this section, we will instead adopt a logarithmic scan over GUT scale parame- 
ters, which favors lower energy soft SUSY breaking parameters over high energy ones. 
Specifically, for a dimensionful parameter a ranging up to a max , and a random number 
x which is sampled uniformly between and 1, we will generate values of a according to 
a = exp(xlog(a ma:r )). 

Applying the log sampling to all dimensionful parameters listed in Sec. |2| (but with 
a linear scan over tan/3), and with the same parameter maxima, we re-plot our results 
in Fig. [|. Here, we see that models with lower values of m^o are sampled much more 
prominently than high values of m^o. In the figure, we see that while bino-like models still 
predict a rather high relic density, now there are some bino-like LSP points with relatively 
small values of m^o and Q^oh 2 ~ 0.1. Meanwhile, the higgsino-like and wino- like bands 



still persist at the same location as in Fig. |T|, albeit with a greater population of points at 
low m^o. We also obtain some mixed dark matter solutions at low m^o. 

To compare against the previous distribution in models (Fig. 0), we again project the 
models as a histogram on the Q^o/i 2 axis in Fig. |I]a). While the minimum in the probability 
distribution still lies nearly at the measured CDM abundance, the minimum has filled in 
somewhat, especially with a mixture of bino-like and higgsino-like solutions. 

For comparison with Fig. Ep)., we apply the requirement of m^o < 500 GeV and gain 
the results shown in Fig. |4p). In this case, the minimum probability has migrated to 
slightly lower values, as the higgsino-like and wino-like solutions with Q^oh 2 ~ 0.1 require 

m^o ~ 800 GeV, and have been rejected. We still have a few bino-like models persisting 
with il^o/i 2 ~ 0.1, although the probability is again quite small (but not as small as in the 
linear sampling case). 



3.3 Comparison of SUGRA-19 results to the mSUGRA model 

In this section, for comparison purposes, we present similar results if we restrict our scan 
to the mSUGRA model with the well-known parameter space 

m , m 1 / 2 , A , tan/?, sign(fi), (3.1) 

with m t = 173.1 GeV as before. We sample m : - 6 TeV, m 1/2 : - 2 TeV, A : -3.5 -»■ 
3.5 TeV and tan/? : 2 — > 60. This expands upon the treatment presented in Ref. plj |, since 
it provides a general scan over mSUGRA parameters. We only present results of the linear 
scan for this case. The results of our scan are shown in Fig. |5[ in the O^-o/i 2 vs. m^o plane. 
We again see that bino-like neutralino solutions populate the high S7,-o/i 2 region, while 
higgsino-like and mixed bino-higgsino dark matter solutions (arising from the focus point 
region) populate the range from Q^o/i 2 ~ 0.01 — 0.3. The nearly pure higgsino solutions 
are almost always lower than the measured DM abundance. The wino-like CDM solutions 
have vanished, since in this case, with gaugino mass unification at Mqut, we always get 
gaugino masses M\ <C M 2 at the weak scale. 

In Fig. ||a)., we project the solutions onto the Q^oh 2 axis. The prevalence of bino-like 
DM solutions is clear, with a maximum around Q^oh 2 ~ 30 — 50. A tail drops off gradually 
at low values of Sl^-o/i 2 , although the tail flattens for a while when it is augmented by the 
appearance of higgsino-like and mixed bino-higgsino DM solutions near the measured value 
of ^o/l 2 ~0.1. 

If we apply the rough electroweak naturalness requirement that m^o < 500 GeV, we 
arrive at the results in Fig. ^6). In this case, many of the higgsino-like solutions with 
S7^o/i 2 ~ 0.1 are now excluded, and the measured relic density lies again at a slight dip 
between the large bino-like maximum and the softer mixed bino-higgsino maximum around 
O^o/i 2 ~ 0.03. By comparing Fig. ||6). for the mSUGRA model against Fig. |2p). for 
SUGRA-19, we conclude that it is actually easier to generate the measured abundance of 
CDM in mSUGRA than in the more expansive SUGRA-19 parameter space. 



4. Scan over PQMSSM with mixed axion/axino cold dark matter 

In this section, we assume the MSSM augmented by the Peccei-Quinn mechanism (PQMSSM) 
is the correct effective field theory below the scale Tr. Thus, the MSSM is supplemented 
by an axion supermultiplet containing the spin-zero and -R-parity even axion a and saxion 
s fields, and the -R-odd spin-^ axino field a. We assume the axino is the LSP, and one of 



the constituents of dark matter. We do not assume a specific axion model, e.g. KSVZ[33] 



or DFSZ[34], nor do we assume a specific SUSY breaking mechanism to give the axino 
mass. We parametrize the low energy axion/axino physics by adopting as free parameters 
the PQ breaking scale f a /N (where TV" is the model-dependent color anomaly), the axino 
mass m a (for which m a < rn^o) the re- heat temperature Tr (which is related to unknown 
inflaton parameters) and the axion field initial mis-alignment angle 6i. 

We will assume a re-heat temperature Tr bounded by mm[10 10 GeV, T^™], where 

Tdcp ~ 10 11 ( 10 i£ GcV ) GeV, so as to be consistent with results from Sec. |3|, which 
avoided overclosure limits from excessive gravitino production in the early universe p5|j. 



This limit also avoids the RTW bound [Q Tr ~ Td cp for which axinos would be in thermal 
equilibrium in the early universe. 3 4 We will also assume m„ ~ 100 keV; for lighter axino 
masses, the thermally produced axinos (see below) would likely constitute warm, rather 
than cold, dark matter p7]]. Thus, we supplement the SUGRA-19 model by an axion/axino 
supermultiplet along with the additional parameters 

• f a /N : 10 9 -»• 10 16 GeV, 

• m a : 100 keV — > m^o, 

• Tr : max[0M GeV,m^o/20] ->■ min[10 w GeV,T dcp ], 

• 6i : —IT —7- 7T. 

The lower bound on f a /N comes from bounds on red giant cooling and supernova 1987a. 
The upper bound on f a /N we take near the GUT scale. The more common bound f a /N ~ 
10 12 GeV comes from assuming 9i ~ 1. The lower bound on Tr comes from requiring 



> 



Tr ~ 10 MeV so that standard BBN can occur[38], and requiring neutralino production at 
temperatures above neutralino freeze-out. 

4.1 Mixed axion/axino relic density 

In this section, we briefly review dark matter production in the PQMSSM, with an axino 
as LSP. Since we assume a value Tr ~ 10 10 GeV, then axinos should never be in thermal 
equilibrium in the early universe [p6| . However, thermal production of neutralinos (now 
assumed to be the NLSP) should proceed as usual, except that now each neutralino will 



3 In this case, if Tr > Td cp , then ma < .2 keV to avoid overclosure due to excessive axino production pq, 
H]. Such light axinos would constitute hot dark matter. 

4 If Tr, ~ f a , then PQ symmetry is restored after inflation, and when the Universe cools, PQ symmetry 
is re-broken, and various domains within a Hubble volume may have different values of mis-alignment angle 
9i. Our upper limit on Tr excludes this scenario (Scenario 1 of Ref. |44|). 



decay via Xi ~ > orf- Since the Xi lifetime is of order 0.1 sec, it is assumed to be BBN-safe. 
This non-thermal contribution to the axino CDM abundance is given by 

q ntp = jm n oh 2 (41) 

since the axino number density from this source is the same as the thermally produced 
neutralino number abundance. 

While the axinos are assumed never in thermal equilibrium in the early universe, 
nevertheless, they can still be produced thermally via radiation off of particles that are in 
thermal equilibrium, i.e. the quarks, leptons, gluons, SUSY particles etc.. This thermal 
axino production mechanism has been calculated in Ref's [40[, [41] and [H5]. The latter 



calculation includes some interaction terms neglected in previous calculations; they find 

ClT p h 2 * 1.24g*F( 9 s)— I R f 1QllGeV V F(^20 5s 2 lnA (4.2) 

y s \ys) QeV 1Q 4 QeV ^ jj n j \ysj y ° 9s y ' 

where g s is the strong coupling evaluated at Q = Tr (e.g. g s (Tn = 10 6 GeV) = 0.932, as 
given by Isajet 2-loop g s evolution in mSUGRA). 

Finally, relic axions will be produced via the vacuum mis-alignment mechanism. In 
this case, we use[43, [P|] 



where f(0i) is the anharmonicity factor. We adopt a recent parametrization of /(#«) from 
Visinelli and Gondolo[f44| given by 



m) 



In 



7^ 



7/6 

(4.4) 



The function f(9i) — >• 1 as Q% — >• and f(9i) — >• oo as 9i — >■ ±7T. 

The total dark matter abundance in the PQMSSM with an axino LSP is then given 
by the sum of the three constituents: 



Qaah 2 = n£ TP h 2 + ftPV + n a h 2 . (4.5) 

4.2 Results of scan over SUGRA-19 with mixed aa CDM 

Our first results from the linear scan over SUGRA-19 model parameters augmented by 
the four PQMSSM parameters are shown in Fig. f|, where we plot models generated in 
the Qaah 2 vs. f a /N plane. Models with dominant axion CDM are in red, while models 
with dominant axino CDM are in blue. The measured CDM abundance is denoted by the 
dashed green line. The densely populated diagonal band which increases as f a /N increases 
is due to mainly axion CDM with Oi ~ 1. The diffuse dotted region at low f a /N, which 
extends up to very large values of Q a ah 2 ~ 10 9 is due to thermally produced axinos. For 
thermally produced axinos, the axino-matter coupling is inversley proportional to / tt /iV, 
and so the coupling is large at low f a /N values, leading to the large thermal abundance. 



As f a /N increases, the axino-matter coupling is suppressed, while the magnitude of the 
axion vacuum energy increases, leading to the case of mixed axion/axino CDM which is 
mainly axions. 

In Fig. ^, we plot models in the £l a ah 2 vs. Tr plane. In this case, we see the models 
are spread over a wide range of Tr values. Large values of Tr give rise to a large rate for 
thermally produced axinos, and so £l a a increases somewhat with increasing Tr. In fact, 
it becomes difficult, but not impossible, to gain a value of Sl a ah 2 ~ 0.1 while maintaining 
Tr ~ 10 6 GeV, as needed for baryogenesis processes such as non-thermal leptogenesis[|4o||. 
Since the axion abundance doesn't depend on Tr, the low Tr points are dominanted by 
mainly axion CDM models. 

In Fig. |9[ we plot models in the Q a ah 2 vs. m„ plane. Here, we see the models are 
spread rather uniformly over many decades of m„, although the tendency is for ^l a a,h 2 
to increase as m„ increases. This is partly due to the fact that each relic axino is more 
massive. 

To compare against results from our thermal neutralino calculations of Sec. ||, we 
project all the models onto the O^/i 2 axis, and present our results as a histogram in the 
number of models generated. We see from frame a), of Fig. [K] that the models are spread 
over many decades of Q a a,h 2 . The measured abundance- denoted by the green dashed line- 
lies on the rising shoulder of the distribution, which in fact peaks around Q a ah 2 ~ 10 2 . If we 
require in addition m^o < 500 GeV as in Fig. g6). as a rough EW fine-tuning requirement, 
then we obtain the results of frame b). Qualitatively, the two plots are similar, and the 
measured abundance still lies on the increasing shoulder of the distribution, which again 
peaks around Q a ah 2 ~ 100. 

4.3 PQMSSM parameters which are preferred by the measured CDM abun- 
dance 

As a final exercise, it may be instructive to see which values of PQMSSM parameters are 
preferred by the measured abundance of CDM in the universe. To this end, we will require 
&<aah 2 to lie between 0.05 and 0.2: i.e. within a factor of two of the measured abundance. 
Then we can project our random scans on the parameter axis. 



In Fig. 11, we project the PQMSSM models onto the f a /N plane. In this case, we see 
that the measured CDM abundance prefers a value of f a /N ~ (1 — 2) x 10 11 GeV, with 
broad tails extending to higher and lower values. This value of f a /N corresponds to that 
which is needed for PQMSSM models with mainly axion CDM and initial mis-alignment 
angle of order 1-2. 



In Fig. 12, we show PQMSSM models with the measured abundance versus axino 
mass rxia- Here, the value of m^ is spread over many decades, with some favoritism for 
small values of rria ~ MeV range. These again correspond to the models with mainly axion 
CDM. 



In Fig. |13, we show PQMSSM models with the measured abundance versus re-heat 
temperature Tr. The Tr values are spread out over many decades. There is some preference 
for low values of Tr in the GeV range, since these models suppress excess thermal axino 
production while allowing for vacuum mis-alignment production of relic axions, which don't 
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depend on Tr. The high Tr tail does extend up to values of Tr as high as 10 6 — 10 GeV, 
which can allow for baryogenesis via non-thermal leptogenesis. 

5. Summary and conclusions 

In this paper, we have reported on a calculation of the thermal abundance of relic neutrali- 
nos in the SUGRA-19 model with 19 free parameters, where all dimensional parameters 
are stipulated at the GUT scale, as suggested by gauge coupling unification, and SUSY 
GUTspj. We find in this rather general framework that bino-like neutralinos populate the 
high Q^oh 2 ~ 10 1 — 10 3 region, while wino and higgsino-like neutralinos tend to populate 
the VL^oh? ~ 10~ 3 — 10 -2 region. The measured CDM abundance sits at a deep minimum 
between these two favored values, especially if one requires rather light sparticles with 
m^o ~ 500 GeV, as suggested by electroweak fine-tuning arguments. In this region, a 
"well-tempered neutralino" is needed; however, these WTNs with the correct relic abun- 
dance are very difficult to generate using GUT scale SUSY breaking parameters. In this 
case, we conclude that it would be a near miracle if SUSY neutralinos were to constitute 
the measured CDM in the universe. This is the opposite conclusion to that suggested by 
the "WIMP miracle", where it is asserted that WIMPs have exactly the right properties 
to constitute thermal relics of the Big Bang. In the case of SUSY WIMPs, this does not 
appear to be so. 

Our results were presented in the context of a linear scan over SUGRA-19 parameters, 
but remain qualitatively valid in the case of a log scan over parameter space. We also 
compare our results with those obtained in the minimal SUGRA or CMSSM model, where 
wino-like x± s don't appear. In the case of mSUGRA, the constrained GUT scale parameter 
choices actually make neutralino CDM somewhat more likely than in the SUGRA-19 case. 

If the strong CP problem is solved by the PQ mechanism in the context of SUSY 
models, then an attractive alternative for dark matter appears: mixed axion/axino CDM. 
A scan over SUGRA-19 parameters, augmented by four extra PQMSSM parameters, shows 
that the mixed axion/axino DM abundance can be generated. While the measured CDM 
abundance doesn't lie at the peak of the predicted ad abundance, it also doesn't lie at 
a minimum. Therefore, in our viewpoint, the mixed axion/axino particles seem a more 
plausible candidate for CDM. 

In the end, the issue will have to be resolved by experiment. A thorough search for 
relic axions with f a /N ~ 10 9 — 10 16 GeV is needed. Currently, only a tiny portion of 
QCD axion parameter space has been explored[46|. A major consequence of our analysis 



then is to motivate our experimental colleagues to consider the possibility of new, more 
encompassing probes of axion dark matter. These should probe much more broadly in 
axion mass, and much more deeply in axion coupling. 
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Figure 7: SUGRA-19 parameter space points from a linear scan where we assume mixed ax- 
ion/axino CDM. There are 4217 points in the figure. 
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Figure 11: SUGRA-19 parameter space points from a linear scan where we assume mixed ax- 
ion/axino CDM. 
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